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<N : Abstract 

I We calculate the contributions from QED processes involving two virtual photons to the cross 

I sections for e^e~ annihilation into two charmonium states with the same C-parity. Generically, 

' the cross sections are three orders of magnitude smaller than those for charmonia with opposite 
C-parity because they are suppressed by a factor of a^/a^. However, if both charmonia have 
quantum numbers J^'~" = 1 , then there is a contribution to the cross section that involves the 

r~| ! fragmentation of each photon into a charmonium. The fragmentation contribution is enhanced by 

^ I powers of -Ebeam/^^o the ratio of the beam energy to the charm-quark mass, and this enhancement 

' can compensate for the suppression factor that is associated with the coupling constants. In 

■ particular, the predicted cross section for J/t/j + J/t/j at the B factories is larger than that for 



PACS numbers: 12.38.-t, 12.38.Bx, 13.20.Gd, 14.40. Gx 



1 



I. INTRODUCTION 



Recently, the Belle Collaboration observed e^e" annihilation into two charmonium states 
at a center-of-mass energy a/s = 10.6 GeV by studying the recoil momentum spectrum of 
the J/'0 [1]. The collaboration measured the cross section for J/z/; + r^c and also found 
evidence for production of J /ip + Xco and J/ip + ric{2S). With more precise measurements 
and additional data, it may be possible to discover the D-wave charmonium states ?7c2(l-D) 
and ip2{^D) at the B factories by studying the momentum spectra of the J/ip and other 
charmonium states. 

The exclusive cross sections for double-charmonium states provide stringent tests of the 
NRQCD factorization method for calculating the cross sections for heavy- quarkonium pro- 
duction [2]. The predictions of the NRQCD factorization method for exclusive double- 
charmonium production reduce essentially to those of the color-singlet model for quarko- 
nium production [3], up to relativistic corrections that are second order in the relative 
velocity v of the quark and antiquark in the charmonium. The nonperturbative factors in 
the NRQCD factorization formula for the cross section in the nonrelativistic limit can be 
determined phenomenologically from electromagnetic annihilation decay rates. Cross sec- 
tions for double-charmonium production can therefore be predicted, up to corrections that 
are suppressed by powers of v^, without any unknown phenomenological factors. 

Recent calculations of the cross section for J/ip + r]c have given results that are about an 
order of magnitude smaller than the Belle measurement [4, 5]. The relativistic corrections 
are large, and they may account for part of the discrepancy [4]. It has also been suggested 
that there may be large nonperturbative corrections to double-charmonium cross sections 
at the 5- factory energy [5]. 

In this paper, we calculate the cross sections for e+e^ annihilation into two charmonium 
states that have the same charge- conjugation parity (C-parity), such as J/ip + J/ip. These 
processes proceed, at leading order in the QCD coupling ccg, through QED diagrams with 
that contain two virtual photons. One might expect that these cross sections would be much 
smaller than those for charmonia with opposite C-parity because they are suppressed by a 
factor of a^/aj.. However, if both charmonia have quantum numbers J^'" = 1 , then there 
is a contribution to the cross section in which each photon fragments into a charmonium [6] . 
The fragmentation contribution is enhanced by powers of £'beam/"^c, where £'beam is the 
beam energy and rric is the charm-quark mass [6]. This enhancement can compensate for 
the suppression factor that is associated with the coupling constants. In particular, the 
predicted cross section for J/t/j + J/ip at the B factories is larger than that for J/ip + rjc- 



II. CROSS SECTIONS 

In this section, we calculate the contributions from e~^e~ annihilation through two virtual 
photons to the production cross sections for double-charmonium states Hi + H2- Charge- 
conjugation symmetry requires the two charmonia to be either two C = —1 states or two 
C — +1 states. We express our results in terms of the ratio R[Hi + H2\, which is defined by 

i;[j/.+/fj= -'f;-^+f . (1) 

where a[e~^e~ — >• = ira"^ / (SE^^^^^) and £^beam = \/^/2 is the beam energy in the 

center-of-momentum (CM) frame. We give the angular distributions dR/dx for each of the 
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helicity states of Hi and H2. The angular variable is x — cos 6, where 9 is the angle between 
the e~ and Hi in the e'^e~ CM frame. 

A. Asymptotic behavior 

The 4 QED diagrams that contribute to the process e+e" — > 7*7* cci + cci at order 
are shown in Fig. 1. The diagrams in Figs. 1(a) and 1(b) contribute only for charmonia 

with quantum numbers J^'" = 1 . such as J/ip, iIj{2S), and ipillD). The procedure for 

calculating the matrix element for e+e^ Hi{Pi) + H2{P2) from the matrix element for 

e^e'^ — > c{pi)c{pi) + c{p2)c{p2) is summarized in Ref. [4]. 

When the energy E'beam is much larger than the charm-quark mass rric, the relative sizes 

of the double- charmonium cross sections are governed largely by the number of kinematic 

suppression factors r^, where the variable r is defined by 

beam 

If we set nic = 1.4 GeV and Ebeam = 5.3 GeV, then = 0.28. 

The asymptotic behavior of the ratio R[Hi + H2\ as r — > can be determined from the 
helicity selection rules for exclusive processes in perturbative QCD [7, 8], which apply equally 
well to QED. For generic charmonia, whose quantum numbers are different from those of 
the photon, only the diagrams in Figs. 1(c) and 1(d) contribute. There is a suppression 
factor for each cc pair with small relative momentum in the final state. Thus the ratio 
R[Hi(Xi)+H2{X2)] decreases at least as fast as as r — > 0. This slowest asymptotic decrease 
can occur only if the sum of the helicities of the charmonia vanishes: Ai -|- A2 = 0. For every 
unit of helicity by which this selection rule is violated, there is a further suppression factor 
of r^. Thus the estimate for the ratio R at leading-order in the QED and QCD coupling 
constants is 

R[Hi{Xi) + H2{X2)] - a2(^;2)'+^i+^^(r2)2+l^i+^2|_ (3) 

The factor of t)^+2^ for a charmonium state with orbital angular momentum L comes from 
the NRQCD matrix elements. At leading order, there may, of course, be further suppression 
factors of r^, and the cross section could even vanish. For charmonia with opposite charge 
conjugation, the estimate for R is given by Eq. (3) with replaced by a^. Thus, the cross 
sections for generic charmonia with the same charge conjugation are suppressed compared 
to those for opposite charge conjugation by a factor of a^/a^ ~ 10~^. 

Production cross sections for charmonia with J^^ = 1 are an exception to this generic 
estimate because they also receive contributions from the diagrams in Figs. 1(a) and 1(b), 
which we will refer to as photon-fragmentation diagrams. These contributions are enhanced 
because the virtual- photon propagators are of order 1/m^ instead of order In the 

amplitude, there are also two numerator factors of rrif. instead of i?bcam! which arise from 
the cc electromagnetic currents. Hence, the net enhancement of the squared amplitude is 
{Ehea.ra/'nT'cy, and the contributions to R can be nonzero in the limit r — > 0. As r — > with 
fixed scattering angle 9, the photon-fragmentation contributions to the cross section factor 
into the cross section for e"'"e~ 77 with photon scattering angle 9 and the fragmentation 
probabilities for 7 — > i^i and 7 — > i72- These fragmentation probabilities are nonzero 
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(2) 




(c) (d) 

FIG. 1: QED diagrams for the process e+e" — > 7*7* —>■ cc\ + cci. The upper and lower cc pairs 
evolve into Hi and H2, respectively. 

at order a only for J^^ = 1 states with hehcities satisfying Ai = — A2 = ±1. The 
contribution to the ratio R for J/-^ + J/'^ has the behavior 

R[Jli^{\i) + J/V^(A2)] ~ Ai = -A2 = ±1. (4) 

This estimate applies equally well to J /ip + V'(25') and ip{2S) + il){2S). For J /ip + il)i{lD) 
and ip{2S) +'4)i{lD) ^ there is an additional suppression factor oiv^ from the NRQCD matrix 
elements. For ipi{lD) + ^i(lZ}), the additional suppression factor is v^. 
Using the QCD version of Eq. (3), we obtain the estimate 

R[J/^P + r^c] ~ al{v'f{r'f, (5) 

which holds generally for S'-wave final states with opposite C-parity. The ratio R in Eq. (4) 
is suppressed relative to Eq. (5) by a factor of (a/os)^ ~ 10~^, but the enhancement factor 
that scales as makes the cross sections comparable in magnitude at the energy of a 5 
factory. 
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B. S'-wave + S-wave 



The cross section for rj,. + rjc receives contributions only from the nonfragmentation dia- 
grams in Figs. 1(c) and 1(d). The angular distribution is 

The total ratio R for rjc + i]c is obtained by integrating over x only from to 1 in order to 
avoid double-counting of the identical final-state particles: 

R[Vc + Vc] = Al - i^i^^^. (7) 

45 m° 

Note that the ratio (7) depends on the charm-quark mass explicitly and also through 
the variable r defined in Eq. (2). One of the factors of (1 — r^)^/^ in Eq. (7) is the nonrel- 
ativistic limit of the phase-space factor PcM/-E'beam, where Pcm is the momentum of either 
charmonium in the center-of-momentum frame. It can be expressed as 



-CM _ AV2(.,M^^,M^J 



(8) 



-^beam ^ 

where X{x, y, z) = x"^ + y"^ + — '^{xy + yz + zx) and is the mass of Hi. 

The angular distribution dR/dx for r]c + i]c{2S) is given by an expression identical to 
Eq. (6), except that one of the factors of {Oi)^^ is replaced by (Oi)^^(2s). Since the two 
final-state particles are distinguishable, the total ratio R for r]c + Vci^^) is obtained by 
integrating over x from —1 to 1. Thus, the total ratio R for 1]^ + rjc{2S) is obtained from 
Eq. (7) by replacing one of the factors of (Oi)^^ by {Oi)j^^{2S) and multiplying by a factor 2. 

The cross section for J /ip + J /ip receives contributions from all four diagrams in Fig. 1, 
including the photon-fragmentation diagrams. The angular distribution dR/dx is 

where the entries of F(Ai, A2, x) are 

F(0,0,x) = rV(l-x')[2 + 3r^-r^ + 4x^(1 -r^)]^ (10a) 
F(±l,^l,x) = (l + x^)(l-x^)[6-7r^ + 4r^-r^ + 4xV^(l-r^)]^ (10b) 
F(±l, ±1, x) = rV(l - x^) [3 - 4r^ + 4r^ - + 4xV(l - r^)]^ (10c) 
F(±l, 0, x) = F(0, ±1, x) = (rV2)[(2 - r'^f{;i - 2r'^ + r'^f - x'^(72 - 228r^ 

+359r'' - 292r^ + 130r* - 32r^° + Sr^^) + 4x^(9 - GOr^ + 122r^ 
-112r^ + 58r^ - 14r^° + r^^) + 16a;V(l - r^){6 - llr^ + llr^ 
-2r^) + 64xV(l-r2)2]. (lOd) 

The total ratio R is obtained by summing over the helicities Ai and A2 and integrating over 
X only from to 1. The resulting expression is so complicated that it is not very useful 
to write it down. The ratio dR/dx for J/ip + ip{2S) is given by an expression identical to 
Eq. (9), except that one of the factors of {Oi)j/^ is replaced by (Oi)^(25) and the range of 
X is from —1 to 1. 
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In the limit r ^ with x fixed, only the transverse helicity states Ai = — A2 = ±1 in 
Eq. (10b) contribute. The angular distribution dR/dx summed over helicities reduces in 
this limit to 

f |J/V.(±l) + J/V(=Fl)l » PU*' (11) 

where P^^jj^ is the probability for a photon to fragment into a J/ip [9]: 

Pi^Jli> = e,7ra ' . (12) 

The coefficient of P'^j^jj^ in Eq. (11) is dR/dx for the process e+e" — >• 77. The differential 
ratio (11) is compatible with the asymptotic form (4). However the integrated cross section 
is enhanced by a factor ln(8/r^), relative to Eq. (4). The logarithmic factor arises because, 
as can be seen from Eq. (9), the potential divergence at x = ±1 in the integral of Eq. (11) 
is cut off at 1 ± a; ~ r*^/8. 

Chang, Qiao, and Wang have pointed out that the cross sections for e'^e~ — > J/ip + j 
is also large if -Ebeam ^ iric [10]. The reason is that this process has a fragmentation 
contribution that factors into the short- distance cross section for e'^e~ 77 and a single 
fragmentation probability P^^j/^. The ratio dR/dx for this process therefore has a finite 
limit as r — > 0, and the cross section is sharply peaked near the beam direction. 



C. S-wave + P-wave 

The cross section for J/ip + he receives contributions only from the nonfragmentation 
diagrams in Figs. 1(c) and 1(d). The angular distribution is 

— [J/ilj{Xi) + hc{X2)]^ 24 <^(^i^^2,x)r (1-r ) ^ — , (13) 

where the non-zero entries of G(Ai, A2, x) are 

^(0,0,0;) = 2{l-x^), (14a) 

G(0,±l,x) = r^(l - + 4x^), (14b) 

^(il,^!,^) = 8(1 -x^), (14c) 

G'(±l,0,x) = 4rV(l + x2). (14d) 

The total ratio R is obtained by summing over the helicities Ai and A2 and integrating over 
X from —1 to 1: 

Rm + he] = ^ r\53 + 22r^)(l - r^f^ (O/).,,.^^,^ _ 
45 



c 



The cross section for rjc + XcJ receives contributions only from the nonfragmentation 
diagrams in Figs. 1(c) and 1(d). The angular distribution is 

f fc + Xc.(AJl = ^-^GjMAI - r^r" ^^ih^^, (16) 
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where the non-zero entries of Gj{X2,x) are 



^0(0,^ 

Gi(±l,a; 
G2{0,x 

G2{±1,X 



= (3/4)(l-x^), 

= (9rV32)(l + x^ + 16x^), 

= (3/4)(l-x2), 

= (9rV32)(l + x^). 



(17a) 
(17b) 
(17c) 
(17d) 
(17e) 



The total ratio R is obtained by summing over the hehcities Ai and A2 and integrating over 
X from —1 to 1: 



R[vc + Xcj] = 



27 



where the coefficients Gj are 



Go — 1, 

Gi = (3/10)(8 + 17r2), 
G2 = (l/2)(2 + 3r2). 



(18) 



(19a) 
(19b) 
(19c) 



D. P-wave + P-wave 

The cross section for he + he receives contributions only from the nonfragmentation dia- 
grams in Figs. 1(c) and 1(d). The angular distribution is 



where the functions H{Xi, X2,x) are 

H{Q,0,x) = 2x^(1 -a;2)[4-6r^ + 3r^- 4x^(1 -r=^)]^ (21a) 

H(±l,^l,x) = hi(x) + hi{-x), (21b) 

H(±l, 0, x) = H(0, ±1, x) = /i2(a;) + h2{-x), (21c) 

//(±1, ±1, x) = (rV2)x'(l - x')[l + (5 - 4x')(l - r')]2. (21d) 

The functions hi{x) in Eqs. (21b) and (21c) are defined by 

h{x) = + x)\l - x'^)[{2 - + 4x{l - xy{l - r^)f, (22a) 

h2{x) = (rV8)(l - xf[{l + 2x){2 - r^) - 8x^(1 + a;)(l - r^)]^. (22b) 

The total ratio R is obtained by summing over the helicities Ai and A2 and integrating over 
X from to 1: 

R[he + he] = (2272 - 3880r=^ + 3300/ - 2158r*^ + 781r^) 

3780 



xr^(l-r^)V2 (Qi)^c(Oi)fec _ (23) 



We have not calculated the cross sections for Xcj + Xcj'- They receive contributions only 
from the nonfragmentation diagrams in Figs. 1(c) and 1(d). We therefore expect them to 
be comparable to the cross section for he + he- 
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E. S'-wave + D-wave 



The cross section for rjc + ?7c2 receives contributions only from the nonfragmentation 
diagrams in Figs. 1(c) and 1(d). The angular distribution is 



f he + VAX,)] = ^ /(A„ .) Ai - rY'^ ^^i^ 



c 



where the functions I{X2,x) are 

l{0,x) = a;2(l-x2)[4-8r2 + 5/ + 4a;2(l-r^)(2-r^)]^ (25a) 
/(±l,a;) = ii{x) +ii{-x), (25b) 
I{±2,x) = i2ix)+i2i-x). (25c) 

The functions ii{x) in Eqs. (25b) and (25c) are defined by 

ii(x) = (3rV4)(l-x)^[2-3r^-2xr^- 8x^(1 + x)(l-r^)]^ (26a) 
i^{x) = (3rV4)(l - x^){l - xY[l + (1 - r^){l + 2xff. (26b) 

Then total ratio R is obtained by summing over the helicities A2 and integrating over x 
from —1 to 1: 

RiVc + Vc2\ = i^dbr (1232 - 3048r^ + 4880/ - 3726r^ + 1439r^) 
22680 

xr^(l-r^)V2 {0i)vA0i),.2 _ ^27) 

The cross section for J/ijj + ipi receives contributions from all four diagrams in Fig. 1, 
including the photon fragmentation diagrams. The angular distribution is 

— [J/^(AO + ^^(A^)] - [A{l-r-){l-x-)+r-f ' ^^^^ 

where the functions J(Ai, A2, x) are 

J(0,0,a;) = 2r^x^{l-x^) [84 + mr^ - 158r^ + lOlr^ - 22r^ + 4a;2(l - r^) 

X (34 - 37r2 + 34/ - 2/) - 32x^(1 - r'^f{2 - /)] ^ (29a) 

J(±l,Tl,a;) = 3+.{x)+3+.{-x), (29b) 

J(±l,0,a;) = 3+o{x)+j+o{-x), (29c) 

J(0,±l,a;) = jo+{x)+jo+{-x), (29d) 

J(±l,±l,a;) = u+{x)+U+{-x). (29e) 
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The functions jx\'{x) are defined by 

= (1/4) (1 - - x f [ - 3(2 - r^)(52 + 16r^ - 38r^ + llr^ + 2r^) 
-24a;r^(2 - r^fi^ - r^) + 4.tV^(1 - r^)(40 - 85r^ + 22/ - 6r^) 
+96a;V(l - r^ f - 32a;V(l - r^ f{2 - 3r^)Y , (30a) 

j+o{x) = (rV8)(l-a;)2[64a;V(l-r2)2 + 128a;V(l-r2)^-8xV(l-r2) 

x(22 + Sr^ + 2r^) - 4x^^(1 - r^)(46 - 13r2 + 8r^) - 2a;(156 - 56r^ 

+4r^ + 38r^ - 13r^) - (2 - r^) (156 - 28r2 - 24r^ + 19r^)] ^ (30b) 

jo+{x) = (rV8)(l + x)2[64xV(l-r2)2-l28a;V(l-r2)2-8xV(l-r=^) 

x(l + 26r^ + 2r^) + 4x^^(1 - r^){52 - Slr'^ + 8r^) - 2a;(156 + 28r2 
-164r^ + 143r^ - 34r^) + (2 - r^){156 - AOr^ + 18r^ + r^)^ , (30c) 

j++{x) = (rV4)(l - x'^) [32a;V^(l - r^f - 4xV^(l - r^){57 - 30r2 + 2r^) - 24^^^ 

X (1 - r^) - a;(156 - 196r2 + 284r^ - 137r^ + 22r^) + 6r^{2 - r^Y ■ (30d) 

The total ratio R is obtained by summing over the helicities Ai and A2 and integrating x 
from —1 to 1. 

We have not calculated the cross section for J/ip + iIj2{^D)- It receives contributions only 
from the nonfragmentation diagrams in Figs. 1(c) and 1(d). We therefore expect it to be 
comparable to the cross section for r/c + r/c2(l-D). 



III. PREDICTIONS FOR B FACTORIES 

In this section, we calculate the cross sections for exclusive double-charmonium produc- 
tion in e+e" annihilation at the B factories. We also analyze the errors in the prediction for 
the J/ip + J/ip production cross section. 



A. Cross sections 

The results in Section II were expressed in terms of the ratio R, which is defined in 
Eq. (1). The corresponding cross sections are 

47rQ!^ 

a[H, + H2] = ^R[H^ + H^]. (31) 
6s 

The ratios R depend on the coupling constant a, the charm-quark mass rric, and the NRQCD 
matrix elements {Oi)h- 

For the charm-quark mass TOc, we use the next-to-leading-order pole mass, which can be 
expressed in terms of the running mass m.c{fhc) as 

rUc = rhcimc) ^1 + ■ (32) 

If we take the running mass of the charm quark to be fhc{fhc) = 1.2 ± 0.2 GeV, then the 
pole mass is rric — 1.4 ± 0.2 GeV. 
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TABLE I: The NRQCD matrix elements {Oi)h- H is a charmonium state. These matrix elements 
are computed in Ref. [4] from electromagnetic annihilation decays, using mc = 1.4 GeV. The errors 
are the statistical errors associated with the experimental inputs only. 



H 


J/r 


rir(2S). v{2S) 


/;,(1P). \rj(lP) 




ID) 


(Oi)ii 


0.335 ± 0.021 


GcV'"' ().] 39 ± 0.010 GcV'^ 


0.053 ± 0.009 GoV^ 


0.095 ± 0.015 


G(^V' 



TABLE II: Cross sections in fb for e+e" annihilation at £'beam = 5.3 GeV into double-charmonium 
states H1 + H2 with C = —1. The errors are only those from variations in the pole mass nic = 1.4± 
0.2 GeV. There are additional large errors associated with perturbative and relativistic corrections, 
as discussed in the text. 







V'(25) 


hc{lP) 






6.65 ± 3.02 


5.52 ± 2.50 


(6.1 ± 0.9)xl0-3 


0.80 ± 0.32 


H2S) 




1.15 ± 0.52 


(2.5 ± 0.4)xl0~3 


0.33 ± 0.13 


K{iP) 






(0.19 ± 0.02) X 10-3 





The NRQCD matrix element {Oi)h can be determined phenomenologically from the 
electromagnetic annihilation decay rate of either i7 or of a state that is related to H by the 
heavy-quark spin symmetry. We make use of the values of the matrix elements that were 
obtained from such an analysis in Ref. [4]. That analysis takes into account the next-to- 
leading-order QCD corrections to the annihilation rate in the cases of the S- and P-wave 
states. It obtains the NRQCD matrix elements for the 25* and ID states from the electronic 
decay rates of the ip{2S) and t/'i(lD), under the assumption that the mixing between the 
25" and ID states is negligible. The values of the matrix elements from that analysis are 
summarized in Table I for the case rric = 1.4 GeV. For different values of rric, the values 
in Table I should be multiplied by (mc/1.4 GeV)^"*"^^, where L = 0, 1, and 2 for 5*- waves, 
P-waves, and D-waves, respectively. 

Our predictions for double-charmonium cross sections are given in Table II for C = — 1 
states and in Table III for C = +1 states. The error bars arc those associated with the 
uncertainty in the pole mass rric only. The small error bars for the r]c + f]c cross section 
in Table III are a consequence of the value of rric being fortuitously close to a zero in the 
derivative of the cross section with respect to rUc- Note that the cross sections in Table II are 
in fb, while those in Table III are in units of 10^'^ fb. All of the cross sections, except those 
for which both of the charmonia are in 1 states, arc at least 3 orders of magnitude smaller 
than the J/ip + J/ip cross section. The cross sections for the 1 states are dominated 
by the photon- fragmentation diagrams in Figs. 1(a) and 1(b). For rUc — 1.4 GeV, the 
photon-fragmentation diagrams contribute 114.5% of the J/t/j + J /if} cross section. The 
nonfragmentation diagrams in Figs. 1(c) and 1(d) contribute 0.9%, while the interference 
term contributes —15.4%. 

The cross sections for production of double-charmonium states with opposite C-parity 
were calculated in Refs. [4, 5]. The cross sections at the P-factory energy in Table II 
of Ref. [4] were calculated using the same NRQCD matrix elements as the cross section in 
Table I. In spite of the suppression factor of a^/a^, the predicted cross section for J/ip + J/il) 
is larger than that for J/ip + rjc by about a factor of 1.8. The kinematic enhancement factor 
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TABLE III: Cross sections in units of 10~^ Po for e^e^ annihilation at -Ebeam = 5.3 GeV into 
double-charmonium states Hi + H2 with C = +1. The errors are only those from variations in the 
pole mass rric = 1.4 ± 0.2 GeV. There are additional large errors associated with perturbative and 
relativistic corrections, as discussed in the text. 



H2\Hi 


Vc Vc{2S) 


Xco(lP) 


Xcl(lP) 


Xc2(lP) 


Vc2ilD) 




1.83 ± 0.10 1.52 ± 0.08 


0.34 ± 0.04 


1.31 ± 0.29 


0.48 ± 0.10 


0.18 lb 0.02 


Vc{2S) 


0.31 ± 0.02 


0.14 lb 0.02 


0.54 ± 0.12 


0.20 lb 0.04 


0.07 zb 0.01 



that scales as more than compensates for the suppression factor of a^/a^. Four of the 
powers of come from the enhancement associated with fragmentation processes. The 
other two powers of come from the fact that the final state J/ip{±l)+r)c violates hadron 
helicity conservation by one unit. The angular distributions do-/d\x\ for rric — lA GeV are 
shown in Fig. 2. The normalizations arc such that the areas under the curves are equal to 
the integrated cross sections 6.65 fb and 3.78 fb. At a; = 0, the differential cross section for 
J/ip + J/ip (normalized as in Fig. 2) is smaller than that for J/ip + rjc by about a factor 
0.66. However the differential cross section for J/ilj+ J /ijj is strongly peaked near the beam 
direction a,t x — ±0.994, where it is larger than that for J/t/j + rjc by about a factor 9.3. 
The reason for the sharp peak is evident from the expression (11) for the limiting behavior 
as r — with x fixed. As we have already mentioned, this peak in the angular distribution 
leads to an enhancement factor ln(8/r"^) in the integrated cross section. 

B. Perturbative corrections 

The perturbative corrections to the cross section for J/'^ + J/V^ have not yet been calcu- 
lated. However the perturbative corrections to the dominant photon-fragmentation diagrams 
in Figs. 1(a) and 1(b) arc closely related to the perturbative correction to the electromagnetic 
annihilation decay rate for J/ijj ^ e+e~, which gives a multiplicative factor 

The perturbative correction to the photon-fragmentation terms in the cross section for J/ip+ 
J/ijj is just the square of the expression (33). If we choose the QCD coupling constant to 
be as = 0.25, which corresponds to a renormalization scale 2m,,, then the perturbative 
correction yields a multiplicative factor (0.79)^ = 0.39. The same perturbative correction 
factor applies to the cross sections for J/V' + '0(25') and '4>{2S) + '4>{2S). This perturbative 
correction factor apphes only to the leading contributions to the cross sections in the limit 
r ^ 0. However, since these contributions are dominant, we conclude that the perturbative 
corrections are likely to decrease the cross sections by about a factor of 3. 

C. Relativistic corrections 

The prescription for calculating the leading relativistic correction for 5"- wave charmonium 
production processes was summarized in Ref. [4]. The leading relativistic correction is 
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FIG. 2: Differential cross sections da /d\x\ for e~'~e~ annihilation at -Ebeam = 5.3 GeV into J /ip+J /ip 
and J /iIj + rjc- The areas under the curves are the integrated cross sections 6.65 fb and 3.78 fb. 
There are large errors associated with perturbative and relativistic corrections, as discussed in the 
text. 



conveniently expressed in terms of a quantity {v'^)h that can be defined formally as a ratio 
of matrix elements in NRQCD. It can also be determined phenomenologically, using the 
mass Mh of the charmonium state as input [11]. 

If we keep only the dominant photon-fragmentation contribution to the cross section 
for J/ip + J/ip, then the relativistic correction to the cross section is closely related to 
the relativistic correction to the electromagnetic annihilation decay rate for J/ip ^ e^e~ . 
The relativistic correction to the photon-fragmentation probability P^^jj^ in Eq. (12) is a 
multiplicative factor 



1 - ^ ko')n 



X 



2m^ 



M 



(34) 



The first factor in the expression (34) arises from the expansion of the electromagnetic 
current in powers of v [12], while the second factor arises from the dependence of the rate 
on the charmonium mass, which can be determined from dimensional analysis. The decay 
rate for J/ip ^ e^e~ is multiplied by the same factor, except that ImdMji^ is raised to 
the power 2 instead of 3. The relativistic correction associated with the phase space for 
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J/il) + J/ij) can be taken into account by replacing the factor (1 — r^)^/^ in Eq. (9) by the 
expression for -Pcm /-E^bcam given by Eq. (8). If all of these correction factors are taken into 
account, both in the decay rate for J/V^ — >• e^e~ and in P^^j/^, then the net effect on the 
prediction for the production cross section for J/ip + J /ip is simply to multiply it by 

2mc \^ -PcMZ-^beam /^i-v 

X (35) 



Mj/^J (l-r2)V2 

For rric = 1-4 GeV, the relativistic corrections decrease the cross sections for J/ip + J/ip, 
J/ip + 1^(23), and ip{2S) + ip{2S) by multiplicative factors of about 0.78, 0.62, and 0.49, 
respectively. The relativistic corrections also decrease the sensitivity to the value of rric- 
The error associated with variations of rric in the cross section for J/t/j + J/t/j in Table II is 
about 41% for nic — 1.4 GeV. After we take into account the relativistic correction in (35), 
this error decreases to about 9%. 

The relativistic corrections to the J/ip + J/ip cross section are significantly smaller than 
and have the opposite sign from the relativistic corrections to the J/ip + rjc cross section, 
which are given in Ref. [4]. For rric — 1-4 GeV, the relativistic corrections to the J/t/j+ric cross 
section are estimated to increase the cross section by about a factor 5.5. The large difference 
in the relativistic corrections suggests that there may be large relativistic corrections not 
only to the absolute cross sections for double- charmonium production, but also to the ratios 
of those cross sections. 



D. Phenomenology 

Recently, the Belle Collaboration has measured the cross section for the production of 
J/ip + Tjc by observing a peak in the momentum spectrum of the J/ip that corresponds to 
the 2-body final state J/ip + rjc [1]. The measured cross section is 

a[J/ip + Tjc] X B[> 4] - (331^ ± 9) fb, (36) 

where B[> 4] is the branching fraction for the r]c to decay into at least 4 charged particles. 
Since -B[> 4] < 1, the right side of Eq. (36) is a lower bound on the production cross section 
for J/ip + Tjc- This lower bound is about an order of magnitude larger than the predictions 
of NRQCD in the nonrelativistic limit [4, 5]. The relativistic corrections are large, and 
they may account for part of the discrepancy [4] . There may also be large nonperturbative 
corrections to double-charmonium production at the S-factory energy [5] . However the large 
discrepancy between the predictions and the measurement is still disturbing. In the Belle 
fit to the J/iJj momentum distribution, the full width at half maximum of the rjc peak is 
about 0.11 GeV. Since the mass difference between the J/t/j and rjc is about 0.12 GeV, there 
are probably J/ip + J/ip events that contribute to the J/ip + rjc signal that is observed by 
Belle. If these were taken into account, they would increase the compatibility between the 
NRQCD prediction of the cross section and the Belle measurement. 

The signal for the production of J/ijj+X, with nix near the mass of the rjc, can be resolved 
into the contributions from J/ip + r)c and from J/ip + J/ip. One way to do this is to measure 
the angular distribution of the production products and fit to a linear combination of the 
angular distributions in Fig. 2. This method tests not only the predicted production rates, 
but also the dominance of the photon-fragmentation production mechanism. Alternatively, 
one can measure directly the rate for the production of J/ip + J/ip with subsequent decay 
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of both J/ip^s into lepton pairs. In the analysis of Ref. [1]. the data set for J/i(j + X, with 
mx near the mass of the r/c, contains only 67 ± 13 events. Since the branching fraction 
for J/ip into e"'"e~ or /^"'"/x" is only about 0.12, that data sample may not be large enough 
to allow one to employ this method. However, as the integrated luminosity increases, the 
double-dilepton signal should become observable. 

The Belle Collaboration also saw evidence for J/ip + Xcoi^P) and J/ip + r]c{2S) events. 
The r]c{2S) was recently discovered by the Belle Collaboration at a mass of M^^(2S) = 
3654 ± 6 ± 8 MeV [13]. Since the mass difference between the '0(25') and the ric{2S) is only 
about 0.03 GeV, any signal for J/t/j+ric{2S) in the J/t/j momentum spectrum is probably also 
contaminated by J/ip + i^{'2S) events. A 3-peak fit to the momentum spectrum of the J/ijj 
by the Belle Collaboration gives approximately 67, 39, and 42 events with an accompanying 
Vc, Xco(l-P); oi ric{2S), with an uncertainty of 12-15 events. The predictions in Ref. [4] for 
the relative cross sections for the production of a J/ip with an accompanying rjc, Xcoi^P), or 
r]c{2S) are 1.00, 0.63, and 0.42, respectively. The observed proportion of events is compatible 
with the NRQCD predictions. 

In summary, we have calculated the cross sections for e+e" annihilation through two 
virtual photons into exclusive double-charmonium states. The cross sections are particularly 
large if the two charmonia are both 1 states. In the absence of radiative and relativistic 
corrections, the predicted cross section for the production of J/ip + J/ip at the B factories 
is larger than that for J/ip + richj a factor of about 3.7. The perturbative and relativistic 
corrections for these two processes may be rather different and could significantly change the 
prediction for the ratio of the cross sections. Nevertheless, the inclusion of contributions from 
processes involving two virtual photons in the theoretical prediction for the cross section for 
J/iJj + Tjf. production is likely to decrease the large discrepancy between that prediction and 
the Belle measurement. 
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Erratum: Exclusive double-charmonium production 
from e+e" annihilation into two virtual photons 
[Phys. Rev. D 67, 054023 (2003)] 
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In Ref. m, we calculated the cross sections for exclusive double-charmonium production from e+e~ annihilation 
into two virtual photons. We made an error in the calculation by omitting a relative minus sign between diagrams 
that differ by the interchange of identical fermions in the final state. There should be a relative minus sign between 
the sum of the Feynman diagrams in Figs. 1(a) and 1(b) and the sum of the Feynman diagrams in Figs. 1(c) and 1(d). 
The error from omitting this sign affects the cross section only if both charmonia have quantum numbers J^^ = 1 . 
A similar error was made in the calculations in Ref. 2]. In what follows, we use results for the J/tp + rjc, J/tp + ric{2S), 
and J/ip + Xcoi^P) cross sections in which this error in Ref. 1^ has been corrected. We show corrections in boldface 
type, except in display equations. 

In the expression for the differential rate dR/dx for J/tp + J/ip production, Eqs. (lOa)-(lOd) should be corrected 
as follows: 

F(0,0,a;) r'^x^{l - x^)[2 + 3r'^ - + 4x^(1 ^ r'^)]'^, (10a) 
F(±l,Tl,a;) = {I + x^){l - x^)[6 - 7r^ + Ar^ - + Ax^r^il ~ r^)f , (10b) 
F(±l,±l,a;) = r^x^il- x^)[3- ir^ +4:r'^ +4:X^r^il-r^)]^, (10c) 
i^(±l,0,a;) = F(0,±l,x) = (rV2)[(2-r2)2(3-2r2+r'')2-a;2(72- 228r2 

+359r'* - 292r6 + 130r^ - 32r^" + 3r^^) + 4x^(9 - 60r^ + \22r'^ 
-\\2r^ + 58r* - Wr^" + r^^) + \%x^r^(\ _ r^)(^ _ \\r^ + \\r^ 
-2r^)^^4x^r^<y\-r^f\. (lOd) 

Eq. (11) and the sentence containing that equation should be deleted. 

In the expression for the differential rate dR/dx for J/4> + production, Eq. (30a) should be corrected as follows: 

J(0,0,x) = 2r'^x^ {1- x^)[84: + 124:r^ -158r^ + 101r'^ -22r^ + 4:X^{1~ r^) 

X (34 - 37r2 + 34r* - 2r'^) - 32x^(1 - ry{2 - r^)]\ (30a) 
Eqs. (31a)-(31d) should be corrected as follows: 

j+-ix) = {l/4){l-x^){l-xf[-3i2-r^){52 + 16r^-38r^ + nr^ + 2r^) 
-24a;r^(2 - r^)^{l - r^) + ix'^r^il - r^){40 - 85r^ + 22r* - 6r^) 

+96x^4(1 - r2)2 - 32x^2(1 - ry{2 ~ 3r^)]\ (31a) 
j+o(x) = (rV8)(l-x)2 [64x^2(1 -r2)2 + 128a;V(l-r2)2-8a;V2(l-r2) 
X (22 + 5r^ + 2r^) - 4x^^(1 - r^){4:6 - 13r^ + Sr"*) - 2a;(156 - SBr^ 

+4:r^ + 38r^ - 13r^) - (2 - r^){156 - 28r^ - 24r^ + I9r^)] ^ (31b) 
jo_^(x) = (rV8)(l + a;)2 [64x^2(1 -r2)2-128xV(l-r2)2_ 8x^2(1 -r^) 
X (1 + 26r^ + 2r^) + 4x^^(1 - r'^){b2 - ilr'^ + 8r*) - 2x(156 + 28r'^ 

-IQAr^ + 143r'^ - 34r**) + (2 - r'^){im - AQr"^ + 18r'* + r^)Y , (31c) 
j++(,x) = (r''/4)(l - x^) [32x^2(1 - r'^f ~ 4x^^(1 - r'^){b7 - iGr"^ + 2r'^) - 24x^2 

X (1 - r^) - x(156 - 196r2 + 284r'* - 137r6 + 22r^) + Qr'^{2 - r'^f] ^ (31d) 
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The error in the calculations produced errors in the numerical predictions for some of the cross sections given in 
Table II of Ref. 0. That table is reproduced here, with the corrected values shown in boldface type. The corrections 
decrease the "5'- wave + S'-wave" cross sections by about 24%. The corrections to the "S'-wave + TD-wave" cross 
sections are almost negligible. 



TABLE II: Cross sections in fb for e^e~ annihilation at iSbcam = 5.3 GeV into double-charmonium states H1 + H2 with C — — 1. 
The errors are only those from variations in the pole mass rric = 1.4 ± 0.2 GeV. There are additional large errors associated 
with perturbative and relativistic corrections, as discussed in the text. 



H2 \ Hi 




^(2S) 


K{\P) 


V'i(l^) 


Jli> 


6.65 ± 3.02 


5.52 ± 2.50 


(6.1 ± 0.9)xl0~^ 


0.80 ± 0.32 






1.15 ± 0.52 


(2.5 ± 0.4) X 10"^ 


0.33 ± 0.13 


hc{lP) 






(0.19 ± 0.02)xl0~^ 





The error in the calculations produced errors in Fig. 2 of Ref. 1]. The corrected Fig. 2 is shown below. The 
normalizations arc such that the areas under the curves are equal to the integrated cross sections 6.65 fb and 3.78 
fb. At a: = 0, the differential cross section for J /ip + J/ip (normalized as in Fig. 2) is smaller than that for J/ip + rjc by 
about a factor 0.66. However, the differential cross section for J/ip + J/ip is strongly peaked near the beam direction 
a,t X = 0.994, where it is larger than that for J/ip + rjc by about a factor 9.3. 
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FIG. 2: Differential cross sections da/d\x\ for e^e~ annihilation into J/iIj + J/i/) and J/ip + rjc at Sbcam = 5.3 GeV. The areas 
under the curves are the integrated cross sections 6.65 fb and 3.78 fb. There are large errors associated with perturbative 
and relativistic corrections, as described in the text. 

In Section III, there are several additional corrections that must be made in the text. 

• For TOc = 1.4 GeV, the photon-fragmentation diagrams contribute 114.5% of the J/iIj + J cross section. The 
nonfragmentation diagrams contribute 0.9%, while the interference term contributes —15.4%. 

• In spite of the suppression factor of a^/a^, the predicted cross section for J/'f/; + J ji^ is larger than that for 
J/f/; + r\c by about a factor of 1.8. 

• The error associated with variations of rric in the cross section for J/ij] + J/'4> in Table II is about 41% for 
rric = 1.4 GeV. After we take into account the relativistic correction in Eq. (36), this error decreases to about 
9%. 

• The predictions for the relative cross sections for the production of a J/ij) with an accompanying r/c, Xco(l^)j 
or ric{'2S) are 1.00, 0.63, and 0.42, respectively. The observed proportion of events is compatible with the 
NRQCD predictions. 
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• The following sentence should be deleted: "If significant fractions of the J/tp + rjc and J/ip + ric{2S) signals 
are actually J/ip + J /ip and J /ip + ?A(2S') events, then the data would be more compatible with the NRQCD 
predictions." 
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